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Oh Context. The luminosity function (LF) is a basic tool in the study of galaxy evolution since it constrains 
O galaxy formation models. The earliest LF estimates in the IR and far-IR spectral ranges seem to suggest 

H 

strong evolution. Deeper samples are needed to confirm these predictions. We have a useful IR dataset, 
which provides a direct link between IRAS and ISO surveys, and the forthcoming deeper Spitzer Space 



Telescope and Akari cosmological surveys, to address this issue. 

Aims. This allows us to derive the 60 /im local LF to sensitivity levels 10 times deeper than before, to 
investigate evolutionary effects up to a redshift of 0.37, and, using the 60/15yum bi-variate method, the 
poorly known 15 fim local LF of galaxies. 

Methods. We exploited our ISOCAM observations of the IRAS Deep Survey (IDS) fields (Hacking & 
Houck 1987), to correct the 60yum fluxes for confusion effects and observational biases. We find indications 
of a significant incompleteness of the IDS sample, still one of the deepest far-IR selected galaxy samples, 
below =^ 80mJy (Mazzei et al. 2001). We have reliable identifications and spectroscopic redshifts for 100% 
of a complete subsample comprising 56 sources with 5 (60yum) > 80 mjy. 

Results. With our spectroscopic coverage we construct the 60yum LF for a sample complete down to 
80 mJy. This LF extends over three orders of magnitude in luminosity, from 9 up to more than 12 in 
log(L6o/^^G)- Despite the fact that the redshift range of our sample exceeds z = 0.3, the V/V„ra.v test gives 
< Vy V„KM >= 0.51 + 0.06, consistent with a uniform distribution of sources. A more direct test, whereby the 
LF was measured in each of four different redshift intervals, does not point out any signature of evolution. 
On the other hand, the rest-frame 15yum local LF we derive, extends up to log(Li5/Lo) = 12 and predicts 
10 times more sources at log(Li5/Lo) =11 than are seen by Pozzi et al. (2004). 



Key words. Galaxies:evolution - Infrared: galaxies, luminosity function - ISM: dust, extinction 
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1. Introduction 



Understanding how galaxies form and evolve is a key goal of modern physical cosmology. A 
fundamental observable of galaxies is their luminosity function (LF) which has long been used 
to constrain galaxy formation models and to quantify star formation and evolution both in lumi- 
nosity and in density. The IR and far-IR spectral ranges are the best to deepen our knowledge of 
this subject since they trace the star formation that is responsible for galaxy formation. In par- 
ticular several satelhte missions, in the past (IRAS and ISO), and in the present (Spitzer, Akari), 
provided and will provide data which will be complementary for a detailed study of the LFs in 
such spectral domains. 

The earliest IR estimate of the LFs, derived from the IRAS data 



Rowan-Robinson et al. 1987 



Saunders et al. 20001) . indicated strong evolution, so that 



LF increases w ith redshift oc ( 1 -h z)^^^ Moreover deep su r veys at 15^m carri ed out 



with ISO (i.e.. 



Elbaz etal. 1999 



Flores 1999; 



Larietal. 2001 



Metcalfe et al. 2003 ) seem 



to require, indeed, strong evolution of ISfim sources (Lagache et al. (2005) and refer- 



ences therein). Several 



Franceschini et al. 2001 



evolutionary models we re developed to explain these results, (e.g., 



Rowan-Robinson 20011) and fit the IR/submilhmiter source counts 



with different degrees of success. Nevertheless none of them is based on a local LF obtained 
from I5um data, since the only available data until recently came from IRAS llfim photometry 



Rush etal. 1993 



Xu 1998 



Fang et al. 19981) . The first attempt to build up the 15//m LF of a 



NEPR subsample was made by Xu (2000). However he said that it must be considered as a 
preliminary work because: i) the sample of galaxies used is a incomplete sample, ii) there is 
a possible misidentification between the sources in the 60/im redshift survey of Ashby et al. 
(1996) and the 15yum sources in his work (see della Valle et al. (2006), for more details), iii) the 
model used to interpret the data treated all IR galaxies as a single population. 

Our IDS/ISOCAM sample overcomes all these issues. It comprises a complete, 60jt/m se- 
lected sample of 56 galaxies in the North Ecliptic Polar Region (N EPR), a subsample of the 



original 98 IRAS Deep Survey (IDS) fields llHacking & Houck 1987 ). 



The IRAS Deep Survey (IDS) sample was defined by co-adding IRAS scans of the 
North Ecliptic Polar Re gion (NEPR), representing more than 20 hours of integration time 
I Hacking & Houck 1987h . It comprises 98 sources with S{6Q/j.m) > 50mJy over an area of 6.25 
square degrees. 

Mazzei et al. (2 001) exploited ISO CAM observations (range 12-18//m) of 94 IRAS Deep 
Survey (IDS) fields llAussel et al. 200(1) . centered on the nominal positions of IDS sources, to 
correct the 60 fim fluxes for confusion effects and observational biases, finding indications of a 
significant incompleteness of the IDS sample below ^ 80 mJy. In della Valle et al (2006) we 
presented spectroscopic and optical observations of candidate identifications of our ISOCAM 
sources. Combining such observations with those by Ashby et al. (1996), we have reliable identi- 



Send offprint requests to: P. Mazzei 



Mazzei et al.: Multi-wavelength study of the IDS/ISOCAM sample.II 3 

fications and spectroscopic redshifts for 100% of the complete subsample comprising 56 sources 
with S (60//m) > 80 mJy. It is the deepest complete IRAS selected sample available and still one 
of the deepest c omplete far-IR selected samples. For comparison, the IRAS Point Source Catalog 



(hereafter PSC; 



Beichma ne^aL^22^ ). com prises about 250, 000 IR sources with a completeness 



limit of 0.5 Jy at 60//m iSoifer et al. 19871) . The deep ISOP HOT surveys, FIR BACK at 170//m, 
(Puget et al. 1999; Dole et al. 2001), and ELAIS at 90//m Joiiver et al. 20od) are all complete 
down to about 100 mJy. Moreover, the 70;um Spitzer catalog of the 8.75 sq. deg. Bootes field 
is flux limited to 80 mJy (Dole et al. 2004) and the Spitzer extragalactic "main" First Look 
Survey, covering about 4 sq. deg., is complete to about 20 mJy at 70 yum (cf. Fig. 2 of Frayer et 
al. 2006), but redshift measurements are available for a substantial fraction of sources (yet only 
72%) merely for^vo^m > 50mJy. 

Thanks to our ISOCAM and optical/near-IR observations our sample, which provides a direct 
link between IRAS and ISO surveys, and the forthcoming deeper Spitzer Space Telescope and 
Akari cosmological surveys is one of the far-IR selected complete samples with the larger 
spectral coverage. In addition to the ISOCAM and to the 60/jm fluxes, most of them ^70% have 
lOOjU m fluxes, the remaini ng 30% with upper limits, and several (^40%) have 25;um fluxes from 
IRAS (iMazzei et al. 20011) . Optical imaging has been akeady performed for 62.5% out of such a 
sample in at least one band, B or R, and for 34% in both the bands, moreover Two Micron All 
Sky Survey (2MASS) data are available for 68% of the sample a nd VLA observations for a large 
fraction of these sources are also available (IHacking et al. 19891) . 

In this paper, which is the second step of our multi-wavelength approach devoted to study 
the evolution of a far-IR selected sample of galaxies on which numerous studies of the far-IR 
evolution of galaxies still rely, we derive the 60//m luminosity function (LF) of such a sample. 
Our sample, ten times deeper in flux density than the PSC catalog, thus less liable to the eff'ect of 
local density inhomogeneity, allows us to investiga te evolutionary effbcts up to a redshift of 0.37, 



Saunders et al. 2000l) . Moreover, we use the 



five times deeper than the PSCz catalog (z< 0.07, 
bi-variate method to translate the 60jum LF to the poorly known 15//m LF. We will compare our 
results with the recent determination of the 15;um local LF obtained by P ozzi et al. (2004) u sing 
the available data on the southern fields, SI and S2, of the ELAIS survey dOliver et al. 200C ). 

The plan of the paper is the following: Section 2 focuses on far-IR properties of our complete 
sample. Section 3 shows our derived 60;um LF, Section 4 presents the 15 //m LF based on the 
bi-variate 60/15yL(m method. In Section 5 there are our conclusions. 
Here and in the following we adopt: A=0.7, Q.h-0.3, Hq-1Q km/s/Mpc. 



' New deep observations of the NEPR are planned with the Akari space mission, also known as the 
InfraRed Imaging Surveyor (IRIS). It will map the entire sky in four far-IR bands, from 50 to 200/jm, and 
two mid-IR bands, at 9 and 20/jm, with far-IR angular resolutions of 25-45 arcsec, reaching a detection limit 



of 44 mJy (5cr sensitivity) in the 50-75yum band jPearson et al. 2004ft . With the Akari orbit, the integration 



time on the NEPR will be particularly high, and, correspondingly, the detection limit significantly deeper 
than average. 
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2. The far-IR properties 



Our complete sample comprises 56 IDS/ISOCAM sources. 25% out of these are beyond z = 0.1, 
12.5% be yond z - 0.15 and o nly 5.3% at z > 0.2 (Paper I). Such a sample is deeper than previous 



estimates llAshbv et al. 1996h . showing a tail extending up to z=0.375, almost 4 Gyr in look-back 
time. 



Our morphological analysis (IBettoni et al. 20061) shows that, although 16% of our sources are 



multiple systems, unperturbed disk galaxies dominate the IDS/ISOCAM sample. One ULIRG, 
3-53A, and two broad Ha emission line galaxies with AGN optical properties (i.e., 3-70A and 
3-96A, see Bettoni et al. (2006) for more details), are also included in the complete subsample. 

Fig.[2(left panel) shows the distribution of the far-IR luminosity (Lf/R, from 42.5 to 122.5 
fim) of our sample where Lp,R = 4; rZ)?(FIR), FIR=1.26xlO"''*(2.58/6o + /loo) W/m^, and /eo 
and /loo are in Jy (iHelou et al. 19881) . In such a figure, as in the following ones, K-corrections 
were derived from evolutionary population synthesis models taking into account dust eff'ects 
I Mazzei et al. 19951) . luminosities were in units of solar bolometric luminosity, Lq - 3.83 x 



10 erg/s. More over, upper limits to flux densities were accounted for by exploiting the Kaplan- 
Meier es timator ( Kaplan & Meier 19581) . Calculations were carried out using the ASURV v 1.2 
package ( Isobe & Feilelson 199oi which implements methods presented in Feigelson & Nelson 
(1985) and in Isobe et al. (1986). The Kaplan-Meier estimator is a non-parametric, maximum- 
likelihood-type estimator of the "true" distribution function (i.e., with all quantities properly 
measured, and no upper limits). The "survivor" function, giving the estimated proportion of 
objects with upper limits falling in each bin, does not produce, in general, integer numbers, but 
is normalized to the total number This is why non-integer numbers of objects appear in the 
histograms of our figures. 

The far-IR luminosity of the IDS/ISOCAM sample extends over 3 orders of magnitude (Fig. 
Q] left panel) with a mean value, log(Lf/R) = 10.2, slightly lower than the mode of the distribu- 
tion, logiLfjn) » 10.5. This value is almost the same as that of the Revised IRAS 6Qfim Bright 
Galaxy Sample dSanders et al. 20031) . and of a normal spiral galaxy, like the Milky Way (Mazzei 
et al. 1992, and references therein). The ULIRG galaxy, 3-53A, emits the maximum far-IR lu- 
minosity of the sample, nearly 100 times higher than the median value. In the same figure (Fig. 
Q] right panel) we present the rest-frame distribution of the /eo/Zioo flux density ratios for our 
sample. This ratio is a measure of the dust temperatu re which gives information on th e relative 



fraction of IR light from new and old star populations (IHelou 1986HMazzei et al. 1992h . Its mean 



value, -0.3, corresponds to a gr ain temperature of about 36 K, consiste nt with the value observed 



for the bulk of IRAS galaxies dSanders et al. 2003 



Soifer et al. 1987h . The Leo/^^ioo ratio cor- 



relates with far-IR luminosity, as expected for the flux density ratio, /eo/Zioo but avoiding any 
redshift dependence, with the most luminous IRAS sources having the largest values of such a 
ratio (see Sanders & Mirabel (1996), and references therein). From the distribution of the rest- 
frame luminosity ratio, Leo/i^ioo we derive a mean value of -0.08, shown as a (red) continuous 
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Fig. 1. Left: The rest-frame far-IR luminosity distribution of the IDS/ISOCAM sample of 56 
galaxies. Here as in Figures|2l|3lS0|9land^2 upper limits are taken into account by exploiting 
the Kaplan-Mayer estimator and by accounting for K-corrections and dust emission using evo- 
lutionary population synthesis models (see text). Right: The rest-frame distribution of the IRAS 
flux density ratio feo/fioOj as revised in Mazzei et al (2001), for the same sample. 



line in Fig.|2] We divide the sample into two subsamples separated about this mean, i.e., a warm 
subsample of 24 IDS/ISOCAM sources with log(Lf,o/Lino) > -0.08, and a co ld subsample of 32 
objects having log(L6o/i'ioo) < -0.08 (iRowan-Robinson & Crawford 19891) . Fig. |3 shows the 
far-IR luminosity distributions of warm and cold sources in our sample. Warm systems entail the 
overall luminosity range with a mode (mean) value, 10.45 (10.3), about 3 (1.5) times higher than 
that of cold sources. 

We use the far-IR luminosity to quantify the star formation rate (SFR). According to 
Chapman et al. (2000), Lj/r = A: x 5 F/? with A: ranging from 1.5 to 4.2, in units of 10^ Lq Mq yr"', 
for Salpeter's IMF with upper and lower mass limit 100 M© and 0.1 Mq respectively and 
Ltir - l.l X Lgo- We adopt k-4.2 by comparing the SFR as defined above with the Hopkins 
et al. (2003) calibration. The more luminous 60 //m sources, i.e., those with higher SFRs, cor- 
respond to the more distant systems, 3-96 and 3-53, two warm sources (Fig.|3}. 

The dust temperatures suggest the presence of two galaxy populations: a spiral population 
with normal dust temperature (i.e., < 36 K), mean SFR ^ 6Mo/yr, mean far-IR luminosity, 
log(Lf/s), ^ 10.1, and mean redshift 0.075, together with a starburst population characterized 
by warm dust temperature (i.e., > 36 K), mean SFR ^ 12M0/yr, mean far-IR luminosity, 
logiLfiR) - 10.3, and mean redshift 0.1 (Fig.|5}. 



3. The 60 fim luminosity function 

The observed flux 5v is related to the rest-frame luminosity by: 



Ly = 4-7TDl Sy lk(L, z) 



(1) 



6 
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Fig. 2. The distribution of the rest-frame luminosity ratio, log L(,qIL\qq, for our 56 IDS/ISOCAM 
sources; red continuous line corresponds to the value log Ldo/Lioo = -0.08 (see text). 




Fig. 3. The rest-frame far-infrared luminosity distribution for the IDS/ISOCAM sample of 56 
galaxies. The separate distributions of warm (24 sources) and cold (32 sources) sub-populations 
are indicated with the (grey) open and (white) hatched histograms, respectively. 

where is the luminosity distance, computed according our cosmological model (©, and 
k (L, z) is the /T-correction defined as: 

(I + Z) iv(l+z) 



(2) 



For the most distant galaxies (z > 0.2) the correction exceeds 20% of the luminosity. 

Fig|6lshows our derived 60yum differential luminosity function, i.e., the co-moving number 
density of sources per logarithmic luminosity interval, (p(L(,Q,z)A logLeo where Lgo - vLy, cal- 
culated using the Schmidt-Eales 1 jVmax estimator proposed by Schmidt (1968) and improved by 
Felten (1976) and Bales (1993). We assume Poisson errors, as tabulated by Gehrels (1986). Lgo, 
i.e., the rest-frame luminosity, has been K-correc ted as specified in §|2l 



According to the VIVmax test llSchmidt 1968h . sources drawn from a population uniformly 



distributed in space should have a mean value of VIV,nax equal to 0.5. Larger (smaller) values of 
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Fig. 4. Left: The SFR distribution for our sample of 56 sources. Right: The SFR distribution for 
warm, 24, and for cold, 32 sources, (grey) open and (white) hatched histogram respectively. 
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Fig. 5. The redshift distribution for warm, 24, and for cold, 32 sources, (grey) open and (white) 
hatched histogram respectively. 



V/V,„ax are the result of a strong increase (decrease) of the co-moving space density of sources 
with redshift. We find < V/V,„„_^ >- 0.5 1 +0.06, consistent with a uniform distribution of sources. 
The IDS/ISOCAM sample is ten times deeper in flux d ensity than the PSCz catalog and 100 
times deeper than the IRAS 60jum Bright Galaxy Sample dSanders et al. 20031) but no signatures 
of evolution arise from our sample. 

The LF we derive is in good agreement with previous wo rks based on the IRAS PSC 



I Saunders et al. 199C : 



Takeuchi et al. 2003 



Takeuchi et al. 2004 . Takeuchi et al. (2003, 2004) 



revised the wor k by Saunders et al. 



(1990) by enlarging their galaxy sample to 15411 galaxies 



from the PSCz dSaunders et al. 20001) with a flux limit of 600 mJy and a redshift range between 
and 0.07. Their analytic fit, shown as a continuous line in Fig.|6l is based on the same parame- 
terization as Saunders et al. (1990): 



exp 



1 



20-2 



(3) 
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Fig. 6. Left.The rest-frame 60;um LF for our sample of 56 IDS/ISOCAM sources (red filled 
squares) compared with LF of IRAS PSCz galaxies from Takeuchi et al. (2003, 2004; blue con- 
tinuous line) and with LF of Spitzer extragalactic First Look Survey from Frayer et al. (2006) 
(open diamonds). The bin size is A log Leo = 0.4. Right: LF split into different redshift bins: 
asterisks (red) for z < 0.05, triangles (green) for 0.05 < z < 0.010, open circles (blue) for 
0.010 < z < 0.015, and stars (magenta) forz > 0.15. 

Table 1. 60 i-im luminosity function: NfWmax in different redshift bins. 



redshift bin 



0.00< z < 0.05 0.490± 0.04 

0.05<z<0.10 0.471±0.08 

0.10<z<0.15 0.444+0.07 

z>0.15 0.315+0.08 



with parameters: a = 1.23 + 0. 04, L. = (4.34 + 0.86) x 10^/?""^ Lp, cr = 0.724 + 0.01, and 



= (2.60+0.30)xlO-^/!'Mpc-^ (iTakeuchi et al. 2003 



Takeuchi et al. 2004h . Our results extend 



over three orders of magnitude in luminosity, from log(L6o/i'0) - 9 up to more than 12. In the 
range where the samples overlap our findings agree with the recent determination by Frayer et 
al. (2006, diamonds in Fig.|6|l, based on a complete sample of 58 sources with S-jq > 50mJy and 
z < 0.3, drawn from the Spitzer extragalactic First Look Survey. 

Despite the fact that the redshift range exceeds z = 0.3, our LF does not show any evidence of 
evolution. 

A more direct test for evolution has been performed by computing LF in four different redshift 
bins assuming no evolution (see Fig.|6l right panel, and TableH*- The results are fully consistent 
with this assumption. 
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4. The 15jum luminosity function 



Deep surveys at 15Atm carried out using ISO (i.e., 



Elbazetal. 199S: 



Flores 199S: 



Lariet al. 2001 



Metcalfe et al. 20031) seem to require strong evolution of 15yum sources starting from redshift 
0.5 (see Lagache et al. (2005) and references therein). Several evolutionary models were devel- 



oped to explain these results (e.g., 



Franceschini et al. 2001 



Rowan-Robinson 20011) trying to fit 



IR/submillimiter source counts with different degrees of success. Nevertheless none of them is 
based on a local LP obtained from 1 5//m data, since the only available data until few years ago 



came from IRAS 12yum photometry ( iRush et al. 1993 



Xu 1998 



Fang et al. 199 



A first attempt to build up the ISfim LF of a NEPR subsample was made by Xu (2000) us ing 



the bi-variate method to translate the 60yL(m local LF of IRAS g alaxies d Saunder s et al. 1990h to 

3 



15/im. Its sample compris es 64 sources detected both at 60/im (iHackin^ & Houck 1987h and at 
15/im ( lAussel et al. 20001) . with redshifts measured by Ashby et al. (1996). Xu fitted the result 
obtained (see its Table 2 and Fig. 3), with a pure luminosity evolution model L oc (1 + f)*'^- 
Nevertheless he said that it must be considered as a preliminary work as discussed in §Q] 

A recent determination of the 15jum local LF was made by Poz zi et al. (2004) us ing the 
available data on the southern fields, SI and S2, of the ELAIS survey dOhver et al. 20001) . Their 
data sample entails 150 galaxies with redshift z < 0.4, excluding sources classified as AGNs (both 
type one and type two). The 15//m LF was calculated with a parametric maximum likelihood 
method. Pozzi et al. (2004) separate spirals from starbursts using optical/mid-IR ratios, assuming 
that Starbursts are the more mid-IR luminous galaxies with, on average, larger log(Li5 /Lr) ratios. 
They estimate the 15yum local LF in the range 7.8 < logiL/Lo) < 10.6. Their findings are that 
the 15fj.m LF of spirals is consistent with no evolution, (< V/Vmax >= 0.55 + 0.03), but that the 
value < y/Vmax >= 0.64 + 0.03 measured for the starbursts suggests that this population is in fact 
evolving. 



4.1. The bi-variate method 

Since our sample is flux limited at 60/im, but not at 15/vm, we used the bi-variate method to 
calculate the 15//m LF of our sample. It was obtained by a convolution of the 60//m LF with the 
Iog(Li5/L6o) distribution: 



0(1 



ogLis) = J <f)(logL(,o)P(log(Li5/L(,o))d(\ogLa)) 



(4) 



where P(log(Lii/L^)) is the conditional probability function that gives the distribution of 
log(Li5) around the mean 15/im luminosity (logiLi^)} at a given 60/im luminosity log(L6o)- 

The distribution per unit interval of the logarithm of the lumino sity ratio is well described by 
Type I Pearson's curves jPearson 1924ilElderton & Johnson 19691) : 



y^yo(l+x/air' (l-x/oz)'"' 



(5) 



where x - log(Li5/L6o), with -ai < x < 02, and origin at the position of the peak of the 
distribution (mode). The values of the parameters are given in Table |2] together with the mean. 
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Table 2. Coefficients and parameters describing the distribution function, Pearson I curve of our 
sample. 



Case 


yo 


ai 


a2 


nil 


nio 


Mode 


Mean 


cr Skew 


Kurt 


-) 

Xy 


All 


82.8 


12.9 


0.43 


83.3 


2.72 


-0.60 


-0.69 


0.24 0.52 


0.52 






16.5 


1.56 


0.04 


10.2 


0.39 


0.04 











Table 3. Parameters describing the LF of our sample. 



Case 




L./lO^La a cr 


All 


5.32 ± .40 


2.00 + .18 1.241 ±.021 0.748 ±.010 



the standard deviation cr, the skewness (fi^/jj^) and the kurtosis (fi4/fil - 3) of the distribution 
(fij is the i-th moment about the mean). The quality of the fit is quantified by the value of 
per degree of freedom (Xv), given in the last column and computed adopting the Levenberg- 
Marquardt method as implemented in Press et al. (1992); cr^ are the errors on the Pearson's 
parameters. 

To perform these analyses we remove from our sample the sources 3-70 and 3-96 with AGN 
properties (see §|2ji and used the parametric solution of the 60/im LF derived by Takeuchi et al. 
(2003, 2004), which agrees well with our results (see §|3 eq.|5}. 

Table |3]lists the parameters defining our 15 /im local galaxy LF using the same parametric 
form as in eq. suggested by Pozzi et al. (2004). The Pozzi LF extends from 7.8 up to 10.6 
in log(Li5/L0), whereas our convolution extends over more than six orders of magnitude and 
100 times deeper (Fig.|S}. Our results agree with those by Pozzi et al. (2004) up to their lumi- 
nosity limit (10.6), however, beyond such a luminosity the two parametric LFs diverge so that at 
log(Li5/L0)=l 1 we expect 10 times more sources than Pozzi et al. (2004). 

4.2. Spiral and starburst populations 

Our sample includes 22 galaxies with rest-frame warm Loo/Lioo ratios and 32 galaxies with 
cold ratios after excluding AGN (see §0. They define two different populations, starburst and 
spiral galaxies respectively, as far as dust properties are concerned. Their distributions per unit 
interval of the logarithm of the luminosity ratios log(Li5/L6o), shown in Fig. |9l (left panel), are 
fitted well by Pearson's curves whose parameters are in Table |3 Their 15 /^m LF fits, i.e., eq. 
Q, are in Fig. |5] (right panel) and the LF parameters are given in Table |5l We find that both 
populations contribute to the faint end of the LF. Spiral galaxies overcome starbursts by less than 
a factor of two. Such a factor is slightly reducing with luminosity, from 1.8 at log(Li5/L0) = 8 
to 1.3 at log(Li5/LQ - 11.8). Our results differ from those of Pozzi et al.(2004), in particular 
predictions concerning starburst population. However they used optical/mid-lR ratios instead of 
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Fig. 7. Distribution of the logLis/Ljo ratio for 54 IDS/ISOCAM sources (see text). The bin size is 
A(log Lis/Lgo) = 0.1. The fitting function is a Type I Pearson whose parameters, optimized following the 
Levenberg-Marquardt method (see text), are reported in Table |2| dashed lines show the same curves with 
Icr error in the parameters. 




Fig. 8. The 15 fim differential LF (red continuous line) for the IDS/ISOCAM sample (54 sources) 
compared with that of Pozzi et al. (2004) (blue dashed line and open triangles). The error bars 
correspond to one cr error of the Pearson's curve parameters. 



far-IR ratios to disentangle starbursts and spirals, assuming that starbursts are the more mid-IR 
luminous galaxies with, on average, larger log{Lis / Lr) ratios. 

To investigate further this point, we repeated our analysis using the same criterion of Pozzi 
et al. to discriminate starbursts from spirals. We used APS catalog l |http://aps.umn.edu| l R-band 
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Table 4. Coefficients and parameters describing the distribution function, Pearson I curve, of 
starburst and spiral populations. 



Case 


yo 


a I 


02 




1112 


Mode 


Mean 


cr 


Skew 


Kurt 


1 

x; 


Starbursts 


35.7 


4.30 


0.14 


9.11 


0.25 


-0.39 


-0.70 


0.25 


1.99 


2.39 


1.70 




12.3 


1.55 


0.02 


3.56 


0.18 


0.002 












Spirals 


43.2 


1.47 


0.32 


4.43 


0.19 


-0.64 


-0.64 


0.23 


0.05 


-0.84 


0.77 


o-,, 


11.6 


0.58 


0.01 


1.89 


0.11 


0.01 













Table 5. Parameters describing the LF of starburst and spiral populations. 



Case 


0./1O- 


3 Mpc-^ 


L./10**Lo 


a 


tr 


Starbursts 


1.89 ± 


.03 


2.96 ± .03 


1.275 ± .001 


0.740 ±.01 5 


Spirals 


3.36 ± 


.07 


2.21 ± .07 


1.236 ±.006 


0.739 ± .027 




Fig. 9. Left; Distribution of the logLis/Leo ratios for starbursts (black continuous line), 22 
sources, and spirals (red dashed line) 32, sources, selected on the basis of their far-IR colors. 
The bin size is A(log Lis/Lgo) = 0.1 and the parameters of the fitting functions. Type I Pearson 
curves, are in Table0]together with one <j errors. Right: The contribution of such galaxy popula- 
tions to the rest-frame 15/im LF (see Table|5}: starbursts (magenta long-dashed dotted line) and 
spirals (green dot-dashed line). Results are compared with those by Pozzi et al. (2004): starbursts 
(red short-dashed line) and spirals (blue long-dashed fine). 

magnitudes for the 37.5% of our sample galaxies that lacked them (see §^, applying a coiTection 
of -0.75 magnitudes to bring the APS zeropoint into agreement with our own common sources. 

Fig. 1101 shows the rest-frame Lis/Lr ratios vs. L15 luminosity. From a least-square-fitting 
procedure we find: 



Lis/Lr = (0.64 + 0.04) X L15 - (6.67 + 0.43) 



(6) 
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log (Lis Lq) 

Fig. 10. The behavior of the log Lis/Lr ratio for warm sources (24, red triangles), and for cold sources (32, 
blue open circles) of IDS/ISOCAM sample; the continuous line corresponds to the least-square fit (see eq. 

Table 6. Coefficients and parameters describing the distribution functions, Pearson I curves, of 
starburst and spiral populations. 



Case yo ai a2 iiii mi Mode Mean cr Skew Kurt xl 

Starbursts 35.11 11.5 0.62 61.2 3.45 -0.63 -0.74 0.29 0.50 0.56 0.86 

a-p 10.7 2.29 0.12 12.3 0.74 0.11 

Spirals 60.3 10.0 0.33 47.8 0.94 -0.60 -0.67 0.21 0.41 -0.34 1.03 

a-,, 13.0 2.61 0.13 12.6 0.47 0.11 



with a dispersion of 0.20dex. Our fit is steeper, 0.64 instead of 0.5, than thatof Pozzi et al. (2004). 
Following Pozzi et al. (2004), we assume Lis/Lr = -0.4 as the nominal separation between spiral 
and starburst populations. Warm and cold galaxies, selected on the basis of their far-IR colors, 
are mixed in Fig.^| Such a criterion selects different galaxies in both the populations. Now there 
are 30 spirals and 24 starbursts in our sample. The distributions per unit interval of the logarithm 
of their luminosity ratios log{Li^/L(,o), are fitted again by Type I Pearson's curves (Fig.[n left 
panel) whose parameters are in Table|6] Table0shows those of their 15 fim LF fits (Fig.[n right 
panel). 

Both such populations contribute to the faint end of the LF. The space densities of starbursts 
are just slightly greater than those of spiral galaxies. 
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Table 7. Parameters describing the LF of starburst and spiral populations. 



Case 


0./1O- 






a 


tr 




Starbursts 


3.82 ± 


.62 


1.81 + .45 


1.245 ± .048 


0.761 ± 


.028 


Spirals 


2.92 ± 


.88 


2.29 ± .85 


1.234 ±.084 


0.742 ± 


.037 





9 10 11 
log (L,5^./Lo) 

Fig. 11. Left; The distribution of the logLis/Ljo ratio for starbursts (black continuous line), 24 sources, 
and for spirals (red dashed line) 30 sources, selected with the same criterion as Pozzi et al. (2004). The bin 
sizes are A(logLi5/L6o) = 0.12 and 0.1 respectively. The parameters of the fitting functions, Type I Pearson 
curves, are in Table|6| Right: The contribution of such galaxy populations to the rest-frame 15pm LF (see 
Table0: starbursts (magenta long-dashed dotted line) and spirals (green dot-dashed line). Our results are 
compared with those by Pozzi et al. (2004): starbursts (red short-dashed line) and spirals (blue long-dashed 
line). 



5. Conclusions 



Combining our observations (Idella Valle et al. 20061) with those by Ashby et al. (1996) we have 



reliable identifications and spectroscopic redshifts for 100% of the complete f ar-IR selected 
subsample comprising 56 IDS sources with S(,o > 80mJy ( iMazzei et al. 200 ll) . The redshift 
distribution shows a tail extending up to z ^ 0.37, in particular ^ 26% of the sources have 
redshifts z > 0.1. 

To fully exploit the potential of this sample, ten times deeper than the IRAS PSC, thus 
less liable to the effect of local density inhomogeneity, for investigating galaxy evolution, 
we calculate the 60jum LF using the l/V„,a.x method. Current estimates are based on rather 
shallow samples. Ev en though our sample is five times deeper in redshift than the PSCz 



Takeuchi et al. 2003 



SaunderseLaL-200C ) used by Takeuchi e t al. (2003), our LF agree with their determination 



Takeuchi et al. 2004^ and with that by Frayer et al. (2006) based on a 



complete sample of galaxies with > 50 mJy drawn from the Spitzer extmgalactic First 
Look Survey. Despite the fact that our redshift range exceeds z - 0.3, and our 60jum LF extends 
up to logLso = 12, whereas that by Frayer et al.(2006) up to logLgo - H, it does not show 
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any evidence of evolution. The V/V„,ax test gives a value consistent with an uniform distribution 
(V/V,nax - 0.51 + 0.06). Moreover, a more direct test for evolution has been performed by split- 
ting the LF in different redshift bins assuming no evolution. The results are fully consistent with 
this assumption. 

We present the bi-variate ISfim LF, one of the few determinations based on ISO data, by con- 
volving the 60//m LF with the luminosity ratio distribution, Lij/Lgo of our sample. This extends 
to luminosity 100 times higher than before. Our result agrees with the recent determination by 
Pozzi et al. (2004) in the common range of luminosity, i.e., from 7.8 up to 10.6 in log(Li5/L0). 
However, above log(Li5/LQ) = 10.6, the two parametric LFs diverge so that at log(Li5/Lo) =11 
we expect 10 times more sources than Pozzi et al. (2004). 

In order to investigate the role of galaxy populations on such a result, we disentangle star- 
bursts and spirals on the basis both of their far-IR dust temperature, and of their Lis/Lu ratios, 
as assumed by Pozzi et al. (2004). Such criteria select galaxies with different dust properties, as 
we point out in § 14.21 Nevertheless, we find that both galaxy populations contribute to the faint 
end of the rest-frame 15/im LF, even though we adopt the same criterion as Pozzi et al. (2004). 
Moreover, in this case, above log(Li5/L0) - 9 our findings are that starbursts and spirals give 
almost the same contribution to the mid-IR LF whereas Pozzi et al. (2004) predict that spirals 
contribute 5-10 times more than starbursts. 
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